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1. Introduction. 

The main claim of this talk is the following: The eigenvalue distribution associated with a 
circular, planar Wilson loop at infinite N undergoes a phase transition of universal character as the 
loop is dilated. The phase transition occurs when the eigenvalue distribution changes from having 
support on an arc with two distinct endpoints to covering the entire unit circle. The transition per- 
sists in the continuum limit for properly regulated Wilson loops, W . For these operators, the short 
distance - long distance cross-over occurs over a scale range that shrinks to zero as N increases 

2. Universality. 

At infinite N, for a small loop, the tail of the distribution at eigenvalues close to — 1 disappears, 
and there is a gap. For finite N the tail is suppressed as e^^^ > 0. For a large loop there is no gap. 
The transition happens at a critical scale where it is governed by a large N universality class. The 
universality class is the same in two, three and four Euclidean dimensions. The two dimensional 
model is exactly soluble and therefore the universality class is fully understood. 

3. The exactly soluble representative of the large universality class. 

A general formulation of the universality class is in terms of a random multiplicative en- 
semble: Consider n independent SU{N) matrices i.i.d. with a measure djJ.. The Wilson loop W 
of unit-less area t = ne^ is represented by a product of the unitary matrices, in any fixed order, 
and evolves in t by diffusion on SU (N) starting from a point source. This is exactly true in two 
dimensions and the relevance of the associated Durhuus-Olesen [^] universality class to higher 
dimensions is an exact version of the old conjecture of "dimensional reduction". 

4. Average characteristic polynomial. 

The average characteristic polynomial Qn of the fluctuating unitary matrix IV at a unit-less 
loop scale t can be defined in any dimension and captures the essential features of the collective 
behavior of the W eigenvalues: 

QN{z,t) = {det{z-W{t))) (4.1) 

We are interested in the region around eigenvalue -1 which affects QN{z,t) mainly around 
z = — 1. In two dimensions the exact result is obtained as e ^ 0, n at fixed t and is critical at 
t = tc = 4. 

4.1 2Ar(z,f) in two dimensions. 

In two dimensions we can calculate QNiz,t) exactly. Applying the Lee- Yang theorem to the 
result we prove that all the zeros of QN{z,t) lie on the unit circle. Figure [I| shows this graphically. 
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4.2 "Double scaling limit" at large A'^ in two dimensions. 

At infinite A'^ tlie logaritlimic derivative of Qn is a two dimensionai eiectric fieid associated 
witli a static cliarge distribution representing the eigenvaiues of W; for one realization of W each 
eigenvalue is counted as a unit charge located at the value of the eigenvalue. When there is no gap 
in the eigenvalue distribution there is a jump in the field at -1 on the unit circle. When there is a 
gap there is no jump. The jump gets smoothed out at finite N ina. universal manner. 

In terms of the characteristic polynomial one captures the smoothing of singularities by going 
to a special limit, in which two new parameters are scaled by N^' and kept fixed as N goes to 
infinity. The region around z = — 1 and t = tc=4 gets infinitely amplified. We restrict our attention 
to the real z axis. One introduces the variables E, and a: 

z = -e'\ - = \+a (4.2) 

with 

2 1 

y= MA <^ = — 7=Ci (4.3) 

We now take N to infinity keeping the variables ^ and a at fixed finite values. The numerical 
prefactors are chosen to make the hmit of QN{z,t) simple. Up to non-universal factors, QN{z,t) 
approaches a function ^{a,^): 

i;{a,^) = J (4.4) 
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Figure 2: Plot of the logarithmic derivative of C, with respect to ^ for different values of the "area" (a) and 
"position" ((^). One sees the smooth remnant of the singularity at = and its dependence on the area. 



d log C,{a,t,) / dt, is similar to a two dimensional electric field and its discontinuity across a charged 
line is smoothed as the line charge density gets depleted in the limit. Figure ^ shows the smooth 
result. 

The meaning of large A'^ universality is that an evaluation of Qn{z,W {t)) also in three or four 
dimensional SU (N) pure Yang-Mills theory, where t labels a dilation parameter and the shape of 
the loop is held fixed, would asymptote to l^{a,^) with the parameters a and ^ defined from y and 
t — tc usmg 

and A^V4 multiplied by some finite numerical amplitudes dependent on the theory 
and on the shape of the Wilson loop. The essence of the claim is that the asymptotic behavior 
will be describable by only two parameters and that the exponents associated with the factors of N 
relating a,^ to will be exactly Vi^ = 1/2,3/4. 

4.3 Convergence to infinite A'^ in two dimensions. 

To test this claim in simulation we first need an indication for how high N ought to be for the 
universal behavior to set in. We plot two cross-sections through the surface of Figure |^ and check 
how rapidly they are approached as N increases in two dimensions, where we have exact formulae 
also at finite A^. The result is shown in Figures ^ and ^ for a Wilson loop that at infinite would 
have a gap and one that would not. For large loops the limit is approached slowly. 

5. Parameter extraction by nonlinear three parameter fits. 

In view of the slow approach to the asymptotic form, we ask how least-square best fits of it 



4 



Phase transitions in Wilson loops 



H. Neuberger 



Convergence to scaling function with increasing N 
d/dy log ( < det [ exp(y/2) + exp(-y/2) W] >) — W = Small loop 
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Figure 3: Wilson loop with a gap. 



Convergence to scaling function with increasing N 

d/dy log { < det [exp(y/2)+exp(-y/2) W ] > ) — W = Large loop 




Figure 4: Wilson loop without a gap. 
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d=2, 

4/t= 



3 parameter nonlinear fits of <der(log(det))> to scaling ^^nction 
b, ranges adjusted centering on b=l, extrapolation by fits to ajj+a^/N "+a,/N 
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Figure 5: Analysis of synthetic two dimensional data. 



to finite A'^ expressions in two dimensions would approximate the values of the critical scale tc and 
the two amplitudes relating y and t — tcto t, and a, when we assume V; to be fixed at 1/2 and 3/4. 
The parameter corresponding to is denoted by A3 and the numerical coefficients multiplying A^^' 
entering a and t, are denoted by Al and by A 1 respectively. In two dimensions the correct values 
areAl =A2=A3 = 1. 

5.1 Two dimensions. 

Figure || shows the result of doing global least square fits to d log l^{a,^)/d^ of synthetic data 
obtained from the exact solution in two dimensions. As indicated in the figure, one extrapolates to 
infinite N by fitting to a series in 1/VN - this form is a consequence of Vi ,2 = 1/2,3/4. 

This exercise shows that one may get close to the right answer even from values of N that are 
practical, but for accuracy at the level of a fraction of a percent one would need to go to A/^ = 100. 
This value is too high for practical Monte Carlo simulations in three or four dimensions. However, 
~ 50 is practical and sufficient for few percent accuracies. 

5.2 Three dimensions. 

We now turn to three dimensions and carry out test runs at some convenient point. Figure ^ 
shows the results. We choose a case where the critical Wilson loop size occurs at an inverse 't 
Hooft coupling b close to unity, so we can just use definitions identical to those in our analysis of 
the synthetic two dimensional data by identifying the parameter t there with ^ here. This is true 
only to 10 percent, but is enough to give us a feel for the order of magnitude of the non-universal 
amplitudes A 1,A2 in three dimensions. 
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d=3, 3 parameter nonlinear fits of der(log(<det>)) to sq^ing function 
ranges centering on b=1.12, extrapolation fits to a^+ai/N "+a,/N 
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Figure 6: Analysis of three dimensional Monte Carlo data. 



Unlike in two dimensions, the data has statistical errors and this further impacts accuracy. 
Nevertheless, one obtains an indication that universality holds with amplitudes numerically quite 
close to their values in two dimensions. This is accidental to some degree and cannot hold in 
general, as in three dimensions there are Wilson loop of shapes unattainable in two dimensions. 

6. Conclusions. 

We conclude that it should be possible to check our main claim in three dimensions but the 
method of analysis ought to be improved. So far, the numerical results in 3D are in agreement 
with large N universality. Amplitudes relating to simple square Wilson loops in three dimensions 
come out surprisingly close to their two dimensional values; perhaps, this is a reflection of the 
approximate dimensional reduction known to hold since the 80's. 
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